Kondo Effects in Carbon Nanotubes: Prom SU(4) to SU(2) symmetry 
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We study the Kondo effect in a single-electron transistor device realized in a single-wall carbon 
nanotube. The K-K' double orbital degeneracy of a nanotube, which originates from the peculiar 
two-dimensional band structure of graphene, plays the role of a pseudo-spin. Screening of this 
pseudo-spin, together with the real spin, can result in an SU(4) Kondo effect at low temperatures. 
For such an exotic Kondo effect to arise, it is crucial that this orbital quantum number is conserved 
during tunneling. Experimentally, this conservation is not obvious and some mixing in the orbital 
channel may occur. Here we investigate in detail the role of mixing and asymmetry in the tunneling 
coupling and analyze how different Kondo effects, from the SU(4) symmetry to a two-level SU(2) 
symmetry, emerge depending on the mixing and/or asymmetry. We use four different theoretical 
approaches to address both the linear and non-linear conductance for different values of the external 
magnetic field. Our results point out clearly the experimental conditions to observe exclusively 
SU(4) Kondo physics. Although we focus on nanotube quantum dots, our results also apply to 
vertical quantum dots. We also mention that a finite amount of orbital mixing corresponds, in the 
pseudospin language, to having non-coUinear leads with respect to the orbital "magnetization" axis 
which defines the two pseudospin orientations in the nanotube quantum dot. In this sense, some 
of our results are also relevant to the problem of a Kondo quantum dot coupled to non-coUinear 
ferromagnetic leads. 

PACS numbers: 75.20.Hr, 73.63.Fg,72.15.Qm 



I. INTRODUCTION 

The first observations of Kondo effect in semiconduc- 
tor quantum dots (QDs) ^i^i^ have spurred a great deal of 
experimental and theoretical activity during the last few 
years. Since these experimental breakthroughs, remark- 
able achievements have been reported, including the ob- 
servation of the unitary limit,- the singlet-triplet Kondo 
effect)^ Kondo effect in molecular conductors^, and the 
Kondo effect in QDs connected to fcrromagnetio^ and 
superconducting reservoirs,^ just to mention a few. 

Recently, Jarillo-Herrero et al. reported perhaps the 
most sophisticated example, namely the observation of 
an orbital Kondo effect in a carbon nanotube (CNT) 
quantum dot (QD).^ In these experiments it was shown 
that the delocalized electrons of the reservoirs can screen 
both the orbital pseudospin degrees of freedom in the 
CNT QD (the K-Y>.' double orbital degeneracy of the two- 



dimensional band structure of graphene) and the usual 
spin degrees of freedom, resulting in an SU(4) Kondo ef- 
fect at low temperatures. In a recent letter fi£ we showed 
that quantum fluctuations between the orbital and spin 
degrees of freedom may indeed dominate transport at low 
temperatures and lead to this highly symmetric SU(4) 
Kondo effect. More recently, Sakano and Kawakamiii 
have studied, using the Bethe ansatz method at zero 
temperature and the non-crossing approximation at finite 
temperatures, the more general case where the quantum 
numbers of N degenerate orbital levels are conserved, and 
found new interesting features of the SU(2Ar)-symmetric 
Kondo effect. Importantly, this is true provided that 
both the orbital and spin indices are conserved during 
tunneling. This poses an interesting question about the 
nature of the nanotube-lead contact because, in princi- 
ple, there is no special reason why the orbital degrees of 
freedom in the CNT should be conserved during tunnel- 
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ing. 

As mentioned, this orbital pseudospin originates 
from the peculiar electronic structure of the nanotube 
(NT)Ai2ii^ The electronic states of a NT form one- 
dimensional electron and hole sub-bands as a result of 
the quantization of the electron wavenumber k± perpen- 
dicular to the NT axis, which arises when graphene is 
wrapped into a cylinder to create a NT. By symmetry, 
for a given sub-band at kj_ — ko there is a second de- 
generate sub-band at k± = — fco- Semiclassically, this 
orbital degeneracy corresponds to the clockwise (O) or 
counterclockwise (O) symmetry of the wrapping modes. 
A plausible explanation of why this degree of freedom 
is preserved during tunneling could be that the QD is 
likely coupled to NT electrodes (the metal electrodes are 
deposited on top of the NT so maybe the electrons tun- 
neling out of the QD enter the NT section underneath 
the contacts) but this issue clearly deserves a thorough 
microscopic analysis about the nature of the contacts. 
The conservation of the orbital quantum number seems 
more likely in the vertical quantum dots (VQD),^^ where 
the orbital quantum number is the magnetic quantum 
number of the angular momentum. 

Here, we take a different route and, assuming some 
degree of mixing in the orbital channel, ask ourselves 
about the robustness of the SU(4) Kondo effect against 
asymmetry in the couplings and/or mixing. 

The rest of the paper is organized as follows: In Sec- 
tion^ we introduce the relevant model Hamiltonian and 
classify different schemes of the lead-dot coupling. These 
different coupling schemes result in different symme- 
tries and hence affect significantly the underlying Kondo 
physics. These effects are analyzed in the subsequent sec- 
tions. Section UTTI presents the analysis with two renor- 
malization group (RG) approaches. In Sec. II VI two slave- 
boson approaches complement the previous results. Fi- 
nally, Sec. Ivl concludes the paper. 



II. MODEL 

A. Nearly Degenerate Localized Orbitals 

We consider a QD with two (nearly) degenerate local- 
ized orbitals which is coupled to reservoirs. As we men- 
tioned before, we have in mind the experimental setup 
of Ref. P where a highly symmetric Kondo effect was 
demonstrated in a CNT QD. However, our description 
could well apply to vertical quantum dot (VQD),^^ where 
the orbitals correspond to two degenerate Fock-Darwin 
states with different values of the angular momentum 
quantum number. Hereafter we denote this orbital quan- 
tum number by m = 1,2. The dot is then described by 



the Hamiltonian 

m=l,2 cr=T,l 

~t~ ^ ^ Ujmyi'fljYKjTlyyi' fjl , (l^) 

{m.a)^{m' ,a') 

where is the single-particle energy level of the lo- 
calized state with orbital m and spin cr, dl^^ {dma) 
the fermion creation (annihilation) operator of the state, 
'T-mcr — dl^^d„ia the occupatiou number operator, Umm 
(to — 1,2) the intra-orbital Coulomb interaction, and 
Ui2 the inter-orbital Coulomb interaction. The effect of 
the external magnetic field parallel to the symmetry axis 
of the system is to lift the orbital and spin degeneracy 
of the single-particle energy levels. We will denote them 
by Aoib and Az, respectively, so that the single-particle 
energy levels have the form 

Cmcr = eo + Aoi-h{Sm,l - Sm,2) + (Az /2){S^^^ ~ Sa^l) ■ (lb) 

The precise values of the Coulomb interactions Umm' de- 
pend on the details of the system, but should be of the or- 
der of the charging energy Ec = e^/2C with C being the 
total capacitance of the dot. In this work we focus on the 
regime where the system of the localized levels is occupied 
by a single electron (J2ma {TT-ma) ~ 1, quarter filling^^) 
and the Coulomb interaction energy {Umm' ~ Ec) is 
much bigger than other energy scales. In this regime the 
Hamiltonian in Eq. (|la() suffices to describe all relevant 
physics of our concern. 

B. Coupling Schemes 

Kondo physics arises as a result of the interplay be- 
tween strong correlations in the dot and coupling of the 
localized electrons with the itinerant ones in conduction 
bands. Naturally, different Kondo effects are observed 
depending on the way the dot is coupled to the electrodes 
and whether or not the orbital quantum number to is con- 
served. Nevertheless, it turns out highly non-trivial ex- 
perimentally to distinguish those different Kondo effects. 
In subsequent sections we will consider different coupling 
schemes between the dot and the electrodes, show how 
different physics emerges, and propose how to distinguish 
them unambiguously in experiments. 

The two leads a — L and R are treated as non- 
interacting gases of fermions: 

-ffa = ^ ^ ^ eafcp °^akt,a°'ak,ia , (Ic) 
k fi— 1,2 cr 

where fi denotes channels in the leads. Without loss of 
generality, we assume that there are two distinguished 
(groups of) channels n — 1 and 2 in each lead. When the 
leads bears the same symmetry as the dot, this channel 
quantum number fi in the leads is identical to the orbital 
quantum number to in the dot and will be preserved over 
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FIG. 1: (color online) Schematic of a representative meso- 
scopic system in question. In (a) Each of the two leads, L 
and R, has two conduction bands (or "modes"), 1 and 2. 
The model with two leads in (a) is equivalent in equilibrium 
to the model in (b) with only one lead. The operators c^^cr 
(^ = 1,2 and a ='\, |) are related to ULk^a and URk^cr by 
the canonical transformation in Eq. l(4j. The wiggly lines in- 
dicate the inter-orbital Coulomb interaction U12 whereas the 
intra-orbital interaction Umm (m = 1,2) is not shown. 



' ^2 {^2^,4.} 

FIG. 2: (color online) Schematics of the SU(4)-symmetric 
Anderson model. 

H in Eq. (Q) is equivalent to _ff = He + Ht + Hd with 



(3a) 



kficr 



Ht= {yf^,mci^^d^a + h.c.^ , (3b) 



where we have performed the canonical transformation 



O'Lkfj.a + CiR.kfj.a- 



(4) 



and discarded the decoupled term ekiib\^^bkiia ■ 

In the following sections we investigate the physics de- 
scribed by the Hamiltonian in Eq. ||2J) and, in partic- 
ular, clarify the role of index conservation in the sym- 
metry of the underlying Kondo regime at low tempera- 
tures. In order to carry out this analysis, we use four dif- 
ferent approaches: the scaling theory (perturbative RG 
approach), the numerical renormalization group (NRG) 
method, the slave-boson mean- field theory (SBMFT), 
and the non-crossing approximation (NCA). 



the tunneling of electrons from the dot to leads and vice 
versa; see Fig. ^a). Otherwise, the orbital channels be- 
come mixed. The most general situation is described by 
the tunneling Hamiltonian 

Ht=Y {Vakt,,^iaa\f,^„d.„,„ + h.c^ (Id) 

akfia 

and the total Hamiltonian is thus given hy H = Hl + 
Hji + Ht + Hjj . 

For the sake of simplicity, we assume identical elec- 
trodes (eLfcp = £Rkfj,) together with symmetric tunneling 
junctions (VLkf^mcr = VRk^Lraa), and ignore their k- and 
CT-dependence of the tunneling amplitudes. In this way, 
we consider a simplified model with Vakfima = Vfj..m/V^ 
and define the widths 

Fm = Fmm j ^mm' — PO^m^m' (2) 

with Vm = Vm.m, whcrc pq is the density of states (DOS) 
in the reservoirs. Then in equilibrium the Hamiltonian 



III. RENORMALIZATION GROUP 
APPROACHES 

The renormalization group (RG) theory provides a 
convenient and powerful method to study low-energy 
properties of strongly correlated electron systems. Here 
we take two RG approaches, the scaling theorjiiSiiLiS 
and the NRG method. ^^iSSiSiiS^ While the scaling the- 
ory is useful for qualitative understanding of the model, 
a more precise quantitative analysis requires the use of 
more sophisticated methods like the NRG method. This 
method is known to be one of the most accurate and pow- 
erful theoretical tools to study quantum impurity prob- 
lems (see Appendix A). 

A. SU(4) Kondo Effect 

We now turn to the case where tunneling processes 
conserve the orbital quantum number; see Fig.|21 In this 
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case, the Hamiltonian reads 

a=L,Rm— 1,2 ka 



+ Hd (5) 



or [see Eqs. lfSa|) and (fSHjl ] 



km a 



E ^™ (4^ 



d\nn^kr. 



Hd- (6) 



From the RG point of view, starting initially with 
nearly degenerate levels, all the localized levels are rel- 
evant for the spin and orbital fluctuations, and, as we 
shall see below, contribute to the Kondo effect. To inves- 
tigate the low-energy properties of the orbital and spin 
fluctuations of the model, we perform the Schrieffer-Wolf 
(SW) transformation and obtain an effective Kondo-type 
Hamiltonian: 



kma'^kma 



H. 



SU(4) 
cff 



(l + 4s-S)(t^r^ + t''T2') 
-(Ji- J2)(s.S)(i-+T-), (7) 



where 



H, 



SU(4) 



J1 + J2 



[s-S + t-T + 4(s-S)(t-T)] (8) 



and the exchange coupling constants Jm {m = 1,2) are 
given by 



1 



1 



(9) 



We note that the Kondo-type effective Hamiltonian in 
Eq. |7J reduces to the SU(4)-symmetric Kondo model 
when Vi = V2 and eia = €20- ■ In this case, orbitals play 
exactly the same role as spins; the former are not distin- 
guished from the latter. 

Under the RG transformation reducing subsequently 
the conduction band width D by 6D, the Kondo-type 
effective Hamiltonian evolves into a generic form: 



leads 



- AzS' - 2Ao,br 



+ 2Ji(s-S) 

+ 2J2(s • S) Q - f 



-[J4 + 4J3(s-S)] (t+T- +t-T+) 



J5f-T'- . (10) 



The level splitting Aorb and Az remain constant under 
the RG transformation: 



dAr, 



dA2 



d\nD d\nD 



0. 



(11) 



The exchange coupling constants Ji (i — 1 , • • • , 5) are 
initially given by Eq. (|5J) and 
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J3 = Ji = VJ^, J5 =2^Ji + J2) , (12) 



which, under the RG transformation, scale as 
dJi 



-2J2- J3(J3-|- J4), 
-2J|- J3(J3-t- J4), 



(13a) 
(13b) 



po^ln D 

dJ2 
podlnD 

'^f' = -JsiJi + J2 + J5) - \Ja{Ji + J2) , (13c) 
PodlnD 2 



dJi ^ T t T T \ 

— Jj — = ~o'^3(-'l + J2) - J4.J5 , 
PodlnD 2 



djr. 



PodlnD 



-3Jq — Ja 



(13d) 
(13e) 



for D 3> Aorb > Az- For D <C Aorb, it is clear from 
Eq. that the orbital fluctuations are frozen and only 
Ji is relevant, which scales as 



PodlnD 



(14) 



It implies that we recover the single-level Anderson model 
for D <C Aorb- Therefore in the remainder of this section, 
we will focus on the case D 3> Aorb- 

It is convenient to define the reduced variables ji = 
Ji/ Ji (i = 2, - - - , 5) and rewrite the RG equations ifT^ 
as 



dh 
dx 

dh 
dx 

dji 
dx 

dh 
dx 



-.12 



= -J3 



'.U 



-J5 



2j1+J3(J3+.74) 

2 + i3(j3+j4) ' 

jsji + h + J5) + + h) 

2 + j3(j3+j4) 
(3j3/2)(l+j2)+i4j5 



2 + i3(i3+j4) 

2 
3 



(15a) 
(15b) 
(15c) 
(15d) 



2 + j3(j3+j4) 

with x — ln{poJi), while Ji obeys the scaling equation 
1 d{poJi) 



{poJiY dlnD 



= -2 - Jsih + Ji) ■ 



(16) 



The RG equations (|15|l have two fixed points: one de- 
scribing the SU(4) Kondo physics 

h = 33 = ji = = 1 (17) 
and the other describing the usual SU(2) Kondo physics 
h = 33 = ji = js = , (18) 
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FIG. 3: (color online) RG flows for different values of r2/ri 
with Fi fixed. 

both with Ji = CX3 as indicated in Fig. 21 Linearizing the 
RG equations (|15|l around the fixed points, one can easily 
show that both the SU(2) and SU(4) Kondo fixed points 
are stable (there is one marginal parameter at the SU(4) 
fixed point). However, as indicated as a dashed semicircle 
in Fig.|3Jb), the radius of convergence is finite while the 
fixed point itself is located at infinity. This implies that in 
priciple, the SU(4) Kondo fixed point cannot be reached 
for arbitrarily small values of 1 — r2/ri. However, as 
illustrated in Fig. Of a) , in the region of physical interest 
for sufficiently small values of 1 — r2/ri, the scaling be- 
havior is essentially governed by the SU(4) Kondo fixed 
point (see also Fig.l^J. More importantly, for sufficiently 
small values of 1 — r2/ri, the SU(2) fixed point governs 
the physics only at extremely low energies. This suggests 
that the SU(4) Kondo signature can be observed exclu- 
sively at relatively higher energy scales (of the order of 
the Kondo temperature), as in the experiment reported 
recently.^ 

At B|| = and Fi = r2 = Fq, the RG equations lO 
reduce to a single equation 

poctlniJ 



0.8- 




FIG. 4: Comparison of the SU(2) and SU(4) Kondo model. 

Comparing this with the corresponding one in Eq. H14|l 
for the usual single-level Anderson model, we note that 
the Kondo temperature is enhanced exponentially: 

T|^(^)^exp(-l/4poJi) (20) 
with respect to the SU(2) Kondo temperature 

T|^(')^exp(-l/2poJi). (21) 

The perturbative RG analysis discussed above, whose va- 
lidity is guaranteed only for pt^Ji <C 1, turns out to be 
qualitatively correct in a wide region of the parameter 
space and provides a clear interpretation of the model. 
To confirm the perturbative RG analysis and make quan- 
titative analysis, we use the NRG method (described in 
Appendix A), the results of which are summarized in 
Figs.21-|31 There the total spectral density 

[E) , (22) 

(7 mm' 

which provides direct information on the linear conduc- 
tance^, is plotted. 

The spectral density shows a peak near the Fermi en- 
ergy, corresponding to the formation of the SU(4) Kondo 
state (see Fig.QJ. The peak width, which is much broader 
than that for the SU(2) Kondo model (represented by the 
dotted line), demonstrates the exponential enhancement 
of the Kondo temperature mentioned above. Another re- 
markable effect is that the SU(4) Kondo peak shifts away 

from Lu = Ef — and is pinned at a; « T^'''^\ This can 
be understood from the Friedel sum rule,^^ which in this 
case gives 5 — tt/A for the scattering phase shift at Ep. 
Accordingly, the linear conductance at zero temperature 
is given by — 4(e^/ft.) sin^ 5 — 2e^ /h. It is interesting 
to recall that the Friedel sum rule gives the same linear 
conductance also for the two-level SU(2) Kondo model. 
Thus, neither the enhancement of the Kondo tempera- 
ture nor the linear conductance can distinguish between 
the SU(4) and the two-level SU(2) Kondo effects. This 
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FIG. 5: (color online) NRG results of the total spectral 
density Ad{E) for different values of coupling asymmetry 
r2/ri. The parameter values are: eo ~ —0.8D, Fi — O.ID, 
Umm' = SD, and Aorb = Az = 0. 
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FIG. 6: (color online) Model with finite mixing between or- 
bital quantum numbers. 



can be achieved only by studying the influence of a paral- 
lel magnetic field in the nonlinear conductance, as shown 
in Ref. ^ 



B. Effects of Mixing of Orbital Quantum Numbers 

To examine the stability of the SU(4) Kondo phenom- 
ena against orbital mixing, we consider the model (see 
Fig.©: 



kcr kma 
+ {""Irnadma + 4..Cfe„-,,) + Hn , (23) 



kma 



where the indices imply 1 = 2 and 2 = 1 and Vq = Vi i 
V2.2 and Vx = Vi,2 = ¥2,1- 



If we rewrite the Hamiltonian in the form 

kcr 



dlAVo 



Ckr. 



krriu 

VxCkrha} + Hd , 



(24) 



it now becomes clear that, in the pseudospin language, 
a finite amount of orbital mixing corresponds to having 
non-coUinear leads with respect to the orbital "magne- 
tization" axis which defines the pseudospin orientations 
m = 1 and m = 2 in the dot. In other words, each 
confined electron (with defined pseudospin) couples to a 
linear combination of pseudospins and, as a result, be- 
comes rotated in the pseudospin space by the angle de- 
fined by tan0 = Vx/Vq. Note that for the maximal 
mixing Vx = Vq, the tunneling electrons lose completely 
information about their pseudospin orientation. In this 
limit, one recovers the spin Kondo physics [with SU(2) 
symmetry] of a two-level Anderson model [see the next 
subsection and Eq. For zero mixing {Vx = 0), the 

model reduces to the SU(4)-symmetric model of Eq. © 
(with tunneling amplitudes which do not depend on the 
orbital index). 

After the RG transformation of the Anderson-type 
model in Eq. (|23|) until the single-particle energy levels 
are comparable with the conduction band width (when 
the charge fluctuations are suppressed), the SW trans- 
formation gives 



ifoff = 1 - 




SU(2) 
off 

4s- s)(t^ + r^) 

+ 2Jx(t"r"), (25) 



where 



H. 



SU(4) 
eff 



= Jo [s • S + t • T 



4(s-S)(t.T)] 

corresponds to the SU(4) Kondo model and 



(26) 



SU(2) 



2 JqS • S(l + 2i^) (1 + 2r=") + Joif" + T^) . (27) 



the SU(2) Kondo model. The exchange coupling con- 
stants Jo and Jx , respectively, are given by 



1 



1 



jo = |yo|^h^ + ^ , j^ = |y^|M_ + _ 



1 



1 



E^ 



(28) 

One can already grasp an idea about the effects of the 
mixing Jx (i-e.. Fx) of the orbital quantum numbers by 
considering the two limiting cases, Jx — (no mixing) 
and Jx = Jo (maximal mixing), of the effective Hamilto- 
nian H25I) . In the case of no mixing (Jx = 0), the effective 
Hamiltonian H25() reduces to the SU(4)-symmetric Kondo 
model in Eq. H26|) . which has been already discussed in 
the previous section: The Kondo temperature is given 
by Tk Z) exp(— I/4J0). When the mixing is maximal 
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FIG. 7: (color online) RG flows in case that there is a flnite 
amount of mixing of the orbital quantum numbers. 



(Jx — Jq), on the other hand, the effective Hamiltonian 
becomes H^^^^^ in Eq. (|?7|). 

Under the RG procedure, the effective Hamiltonian 
(|25|l transforms to the general form 



H,e = Jis • S + [J3(i"r^) + 4J2(s • S){t'T')] 

+ [Jiif'T'' + tyry) + 4j6(s • s)(i^r^ + f^r^)] 
+ [j^it^'T'' ~ tyry) + AJris ■ s)(t^T^ - f^r^)] 

+ [J9+4J8(s-S)(r+T^)] , (29) 



where the exchange coupling constants are initially given 

by 



Ji — J& — Jo 1 J5 = J7 = Jx 



Under the RG transformation, they scale as 
dJi 



podlnD 

dJ2 
podln D 

d-h 
Pod hi D 

dJ^ 
PodlnD 

dJz 
PodlnD 

dJe 
Podln D 

dJ-r 
PodlnD 

dJs 
PodlnD 



J'2 8Jq , 



— — 2 Ji J2 — 2J4 Jg + 2J5 J7 , 



—J3J4 — 3J2 Jg , 

J3J5 + i'h'h , 

—J2J4 ^ 2 Ji Jg — J3 Jg — 4j| , 
J2J5-2J1J7 + J3J7-4JI, 

— 2 Ji Jg — 2Jg Jg — 2J7 Jg , 



and 



dJg 
dlnD 



0. 



(31a) 
(31b) 
(31c) 
(31d) 
(31e) 
(31f) 

(31g) 
(31h) 

(32) 



As before [see Eq. (|15ll ]. it is convenient to work with the 
reduced coupling constants ji = J^ / Ji . In terms of these 
reduced constants, the RG equations read 



(33a) 
(33b) 
(33c) 
(33d) 
(33e) 
(33f) 
(33g) 



dj2 _ 


-J2^ 


2j2 


+ 2j4 je - 2j5i7 


dx 




+ 2jI + 2j2 + 8jg2 


djs _ 


-J3 ^ 




- jI + 3j| - 3j2 


dx 




+ 2ji + 2j2 + 8j| 


dji _ 


-J4^ 




dx 




+ 2j| + 2j2 + 8jl 


djs _ 


-35 - 


izk + 5j237 


dx 




+ 2jI + 2j2 + 8j| 


dje _ 






- 2^6 + jajg + 4j| 


dx 




+ 2jI + 2j2 + 8j| 


djr _ 




kk - 


- 2^7 + k37 - 4j2 


dx 




+ 2j| + 2j2 + 8jg2 


djs _ 


-J8^ 


2j8 


+ 2jg jg + 2j7j8 


dx 




+ 2jI + 2j2 + 8j| 



together with 

1 d{poJi) 
{poJif dlnD 



-(l+j| + 2j| + 2j72 + 8j| 



(34) 



Jg^Js^ VJoJx ■ (30) 



The RG equations lO again have two fixed points, 
one associated with the SU(2) Kondo effect and the other 
with the SU(4) Kondo effect; see Fig. □ The RG flow 
diagram in Fig. [7|is reminiscent of that in Fig.O Both 
fixed points are stable. However, since the radius of con- 
vergence of the SU(4) Kondo fixed point is finite, the 
SU(4) Kondo fixed point cannot be reachable even for 
arbitrarily small mixing Vx [see Fig. EJb)]. However, in 
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FIG. 8: NRG results for the effects of the finite mixing of 
the orbital quantum number m. The parameter values are: 
eo = ~0.8D, To = 0.08D, Umm' 



2T\ 



51/(4) 



8D, Aorb = 32r^^'''', and 



the region of interest, the physics is essentially governed 
by the SU(4) Kondo fixed point for sufficiently small Vx 
[see Fig. [7fa)]. Therefore, the SU(4) Kondo physics are 
in principle unstable against both the orbital quantum 
number anisotropy 1 — r2/ri and the orbital mixing Tx- 
For sufficiently small values of those, however, the SU(4) 
Kondo physics still determines the transport properties 
except at extremely low energy scales. As addressed al- 
ready, this suggests that to observe indications of the 
SU(4) Kondo physics exclusively, one has to investigate 
the properties at relatively higher energies (of the or- 
der of the Kondo temperature). This is confirmed and 
demonstrated in the NRG results summarized in Fig. |S1 
We will also see below that there is no way to distinguish 
the two-level SU(2) Kondo physics and the SU(4) Kondo 
physics experimentally by means of linear conductance. 



C. Two-Level SU(2) Kondo Effect 



As pointed out in the previous subsection, at the max- 
imum mixing (Vq = Vx) the physics becomes that of 
the two-level SU(2) Anderson model (see Fig.^. In this 
case, the only degree of freedom which is conserved dur- 




{d2a,4a} 

FIG. 9: (color online) Schematic of the two-level SU(2)- 
symmetric Anderson model. 



ing tunneling is the spin and the total Hamiltonian reads 

B. k(j 

XI ^™ ("i/ca^™'^ + dlaa'^ocka) + Hd (35) 



a—L,R k(7 



akmcr 



or equivalently [see Eq. (0J] 

ka 

)+Hd. (36) 

kma 

As the scaling theory of the Kondo-type Hamiltonian 
obtained from the two-level Anderson model has been 
developed in detail in Refs. and|2^ we here focus on 
the first stage, which highlights the difference between 
the two-level SU(2)-symmetric Anderson model and the 
SU(4)-symmetric Anderson model. Finally, the physical 
arguments based on the perturbative RG theory will be 
examined quantitatively by means of the NRG method. 

As we integrate out the electronic states in the ranges 
[-D,-{D-SD)] and [D-8D, D] of the conduction band, 
the dot Hamiltonian (|la|l evolves as 



(37) 



Hd = f^madl^^dma - ^ X] {d\^d2cr + dl^di. 

(m,cr):^(m',cr') 



with other terms in the total Hamiltonian (|35|l kept un- 
changed. Notice here the appearance of a new term in 
t, i.e., a direct transition between the two orbitals m = 1 
and 2. Scaling of the parameters Cma and t are governed 
by the RG equations 



dcr, 



= — r„ 



d\nD TT 



and 



respectively. 



dt 



d\nD TT 



(38) 



(39) 
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D/T 

FIG. 10: (color online) Scaling of the single-particle energy 
level ed, to be compared with emo- in Eq. of the single- 

level Anderson model, = ed{D — F) is a scale-invariant 
quantity. 

The RG equation (|38() for the single-particle energy 
levels emu is the same as that in the usual single-level An- 
derson modeliiiS (the corresponding RG flow diagram is 
shown in Fig. I1(J|) . However, due to the direct transition 
t emerging from the RG equation H39(l , ema are not rele- 
vant to the Kondo effect [they are not the eigenvalues of 
Hd in Eq. (|37|l ]. To find the relevant energy level(s) di- 
rectly involved in the Kondo effect, one may diagonalize 
H]j in Eq. H37|) by means of the canonical transformation 







cos(0/2) 


sin(6'/2)' 










-sin(6'/2) 


cos(6l/2) 




d2a 



(40) 



where the angle (j) is defined by the relation 

tan 6*= — - — . (41) 
£2 ~ ei 

The dot Hamiltonian in Eq. (|37|l then takes the form 

M=± f {m,a-)^{m' .it') 

(42) 

with 

e,=± = l{ei + e2)TlVi^i'-^2r+t'- (43) 

At the same time, the canonical transformation in 
Eq. H4U|) also changes the coupling term in the total 
Hamiltonian in Eq. H36|l to 

= E E {^Ld,. + dl^Ck.) (44) 
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FIG. 11: (color online) Scaling of the two-level SU(2)- 
symmetric Anderson model. The arrowed lines indicate RG 
flows of the effective single-particle energy levels e± [see 
Eq. I|43|l ] and the widths of the shadowed regions around e± 
the RG flow of F± [see Eq. ^]. The mean of Fi and F2 is 
denoted by Fq. 

with V± defined by 



'v+ 




cos(6l/2) sin(6'/2)' 




vr 


V- 




-sin(6'/2) cos(6l/2) 




V2 



(45) 



Accordingly, the tunneling rates r± = 7rpo|V±|^ of the 
effective orbital levels e±^a are given by 

r± = i(ri+r2)± Vrlr^sin0+i(ri-r2)cos0. (46) 

Figure 1111 shows the scaling of e± (arrowed lines) and 
F-i- (widths of the shadowed regions around e±), governed 
by Eqs. gSl, and Note that the effective 

single-particle energy levels e± always repel each other )2i 
and the hybridization r+ (r_) of the lower (upper) level 
e+ (e_) always increases (decreases). Essential in this 
scaling property of the two-level Anderson model is the 
direct transition t between the orbitals m = 1 and 2, 
mediated by the conduction band. 

The scaling of e± and F± stops when the lower level 
£+ becomes comparable with D (e+ ~ D); see Fig. ITTI 
Then the charge fluctuations are highly suppressed and 
the occupation of the lower level becomes close to unity 
((n+) w 1). Therefore, only the lower level e+ gets in- 
volved in the Kondo physics, and hence the resulting 
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Kondo effect is identical to the usual SU(2) Kondo ef- 
fect. To be more specific, let us consider the two limiting 
cases, |ei — £2] S> Fq and |ei — £2! ^ Tq, assuming 



|ri-r2| «ro = (ri + r2)/2. 

Since t ^ Tq, one has 

^_K/2, |£i-£2|«ro; 
~|0, |£l-£2|»ro. 



(47) 



(48) 



or equivalently. 



2ro, |£i-£2|«ro: 
To , |£i - £2! > To . 



(49) 



This implies that when the two orbital levels are nearly 
degenerate (|£i — £2] ^ Fq), the Kondo temperatureiii^ 
is enhanced exponentially: 



1 



7r£o 

2ro 



(50) 



[with £0 = (£1 + £2)72] compared with the single-level 
case (i.e., \ei - £2] > Tq) 



T° ~ -v/rVDexp 



To 



(51) 



In the limit of nearly degenerate levels ( |£i — £2 1 ^ Fq) , 
the upper level £_ is located at distance smaller than r+ 
from the lower level £+ [(£_ — £+) < r+; see the up- 
per panel in Fig. ^] and the transition from £+ to £_ 
is allowed in general. Indeed, this effect can be taken 
into account by a proper SW transformation includ- 
ing both levels and scaling of the resulting Kondo-like 
Hamiltonian;SiS& and gives rise to a bump structure at 
u; — Acff above the Fermi energy Ep of the leads, with 
Aeff given hy^ (with ei^ = £2^ initially) 



A 2Fo 

Acff In — 

TT Fo 



(52) 



in the single-particle excitation spectrum Ad(Lo) in 
Fig. see below. 

Again, all the interpretations made above on the ba- 
sis of the perturbative RG are confirmed with the NRG 
method. Figure El shows the total spectral density 
Ad{E). One can see that as Aorb increases with Az = 0, 
the Kondo peak gets sharper, i.e., the enhancement of 
the Kondo temperature Tk in Eq. H50|) diminishes for 
Aorb > Fq; see Fig. I12r a'). Notice that the bump above 
the Fermi energy originates from the excitation via the 
transition from the lower level £+ to the higher one £_, 
and is thus located aX E = Acff [see Eq. (|52) . When we 
allow Az finite as well, the Kondo peak then splits into 
two because of the Zeeman splittingiSSiSS 
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' 1 ' 
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FIG. 12: (color online) Total single-particle excitation spec- 
trum Ad{u}) with (a) only the orbital degeneracy lifted 
(Aorb 7^ 0, Az = 0) and (b) both the orbital and spin de- 
generacies lifted (Aorb,Az 7^ 0). The short vertical arrows 
indicate the transition from e+ to e_, whose excitation en- 
ergy is given by Aefi [see Eq. 1521 1. The parameter values 
are: eo = -0.8D, To = O.ID, and f/™' = 8D. 



IV. SLAVE-BOSON TREATMENT 



In order to confirm our previous results and obtain 
analytical expressions for intermediate mixing, we also 
use slave boson techniques. In particular, the SBMF 
approach, which provides a good approximation in the 
strong coupling limit T <^ Tk, allows us to obtain an- 
alytical expressions for the Kondo temperature and the 
Kondo peak position for arbitrary mixing. Our SBMFT 
results are complemented with the NCA, which takes into 
account both thermal and charge fluctuations in a self- 
consistent manner. 



At equilibrium it is convenient to change into a rep- 
resentation in terms of the symmetric (even) and anti- 
symmetric (odd) combinations of the localized and de- 
localized orbital channels. Thus the even-odd trans- 
formation consists of afe,i(2),cr = (cfe^o- i ick^a)^/^ and 
di(2)a = {deer i idou) / lu this basis the Hamiltonian 
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in Eq. Q reads 

cr,i/— e,o a,u—e,o 

+ ■^011 (4„,.^o- + h.c) , (53) 

where, again we have taken Vq — V14 = V2,2, Vx — 
Vi,2 = ^2,1, Um,m' = U, and ei,^cr = eo,o-- The occu- 
pation per channel and spin is given by n,ycr = d^^^d^^ 
and the total occupation per channel is Ui, = '^„n^„. 
The tunneling amplitudes for each channel are given by 
Ve = Vq + Vx and Vo = Vq — Vx- In order to nor- 
malize the total tunneling rate, we take, for the diago- 
nal and off-diagonal tunneling amplitudes, Vq = V cos (j) 
and Vx = Vsmcj), namely Ve = V{cos4> + sin0) and 



Vo = V{cos(f)— sincf)). Notice that one needs (j) € [0, 7r/4] 
in order to have always Vo positively defined. There ex- 
ist two very different situations, namely (i) (f) — 0, where 
there are only tunneling processes that conserve the or- 
bital index, and (ii) — 7r/4, where the mixing and no 
mixing tunneling amplitudes are the same. 

Now we write the physical fermionic operator as a 
combination of a pseudofermion operator and a boson 
operator as follows: d^o- — h^fua, where the pseud- 
ofermion operator f^^„ annihilates one "occupied state" 
in the vth localized orbital and the boson operator cre- 
ates an "empty state" . Quite generally, the intra- /inter- 
Coulomb interaction is very large and we can safely take 
the limit U ^ co. This fact enforces the constraint 
X^i/cr flafi^o- + i>^b — 1 that prevents the accommodation 
of two electrons at the same time in either the same or- 
bital or different orbitals. This constraint is treated with 
a Lagrange multiplier. 



h.c. 



e,o 



(54) 



Notice that we have rescaled the tunneling amplitudes 
Ve(o) VNVe{o) according to the spirit of a 1/A^- 
expansion (where N is the total degeneracy of the lo- 
calized orbital). 

Our next task is to solve this Hamiltonian, which is 
rather complicated due to the presence of the three op- 
erators in the tunneling part and the constrain. In order 
to do this we employ two approaches that describe two 
different physical regimes. The first one is the SBMFT 
approach which describes properly the low-temperature 
strong coupling regime. This SBMFT provides a good 
approximation in the deep Kondo limit, namely in case 
that only spin fluctuations are taken into account. The 
NCA, on the other hand, takes into account both ther- 
mal and charge fluctuations in a self-consistent manner. 
It is well known that the NCA fails in describing the low- 
energy strong coupling regime. Nevertheless, the NCA 
has proven to give reliable results at temperatures down 
to a fraction of Tk- 



A. Slave-boson Mean-Field Theory 

We begin with the discussion of the mean-field approx- 
imation of Eq. H54|) . The merit of this approach is its 
simplicity while capturing the main physics in the pure 



Kondo regime. It has been successfully applied to inves- 
tigate the out-of-equilibrium Kondo effecli^iiSSiSi^ and 
double quantum dots,^'''"^^ ''^''^'^'''^ just to mention a few. 
This approach corresponds to taking the lowest order 
^(1) in the 1/N expansion, where the boson operator 
b{t) is replaced by its classical (nonfluctuating) average, 
i.e., b{t)/y/N {b)/^/N = b, thereby neglecting charge 
fluctuations. In the limit N — + oo, this approximation 
becomes exact. The corresponding mean-field Hamilto- 
nian is given by 

e,o i^,a 

+ E + H 

+ A(^^0.,. + iV|6p-l^ ' (55) 

where Vi, = V^b — Vi,{b) and eo,o- — eo.cr + A are the 
renormalized tunneling amplitude and the renormalized 
orbital level, respectively. The two mean-field parameters 
b and A are to be determined through the mean-field 
equations, which are the constraint 

YifUt)^At))+N\b\^ = l (56) 
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and the equation of motion for the boson field 



-X ~ Eq, these equations obtain the forms: 



vAc'At)UAt)) + mbV^Q- (57) 



The Green function for the v {= e,o) locahzed orbital 
and the corresponding lesser lead-orbital Green func- 
tion are given by G<„{t - t') = -i{flAt')f^,At)) 

and G^^j^^^t -t') = -i{cl^^^it')fu,ait)) , respectively. 
Expressing the mean-field equations in terms of these 
noncquilibrium Green functions, we obtain Eqs. (|56|l and 
()57|1 in the form 



a 



0. 



(58a) 
(58b) 



In order to solve the set of mean- field equations, we 
proceed as follows: First, we employ analytic continua- 
tion rules to the equation of motion for the time-ordered 
Green functions Gl^t-t') = -i{Tc{fUt')fUt)}) and 

Gl,.;k.Jt-t') = -^{Tc{cl jt')UAtm, where Tc de- 
notes the time-ordering operator along a complex time 
contour. Second, we use the equation of motion tech- 
nique to relate the lead-orbital Green function with the 
orbital Green function. Finally, we rewrite the mean-field 
equations in the frequency domain (taking £0,0- — So)'- 



Mb? 



-E 



de 
2ni 



G<,(e)(eo-eo) = 0. 



(59a) 
(59b) 



The integrals in Eq. H59|l can be carried out analytically 
by introducing a Lorentzian cutoff p{e) = D{e^ + D^)~^ 
for the DOS in the leads and the lesser orbital Green 
function G<(e) = 2if^/(e)/[(e - en)^ + f^] with f^, = 
Fi/j&p and the Fermi distribution function /(e): 





[ln(^^^° + ^^^)| 


2 ^ 
H Im 


[ln(^^~" + ^^°)| 


ttN 




ttN 





2r, 



1 

N ~ 
-Re 



\b? 



(60a) 



ttN 
= -X\b\ 



[ln(^^~° + ^^^)| 


^ 77 Re 


[ln(^^~° + ^^°)| 




ttN 





(60b) 



Im 



In 

TT 

2 ' 



£0 + i^e 

D 



Im 



In 



£0 + iTp 
D 



FeRe 



In 



£0 + iFe 
D 



FoRe 



In 



£0 + iTp 



Treo 
2 ' 



(61a) 



(61b) 



where Fi, = T^/N is the original rate for the v (= e, o) 
channel. Using the parametrization Ve — ^(cos 0-|-sin </>) 
and Vo — V{cos(j) — sin0), the tunneling rates read Fg = 
ttV'^p{1 + sin 2(/i) F(1 + sin 20) and To = p = nV'^p{l - 
sin 20) = F(l - sin 20). Taking sin 20 = /3 with P e 
[0,1] (notice that < sin 20 < 1 for e [0,7r/4]), we 
parametrize the even and odd rates as Fe = {1 + P)T and 
Fo = (1 — /3)F, respectively. Accordingly, the case /? ^ 
accounts for the process where even and odd channels are 
not coupled equally to the lead electrons or equivalently, 
the process where the orbital index is not conserved. In 
terms of the new notation, the mean-field equations can 
be written in a compact way: 



, (62) 



In 


£0 + «f e 


+ ln 


£0 + «f 


+ In 






D 




D 







TT TTEq 

i 1 

2 Fe + Fo 



or equivalently, 
[ea+iTe][eo+iTo] 



£^ + f? 



-0 ^ o 



(63) 



the real and imaginary parts of which read 



13/2 



= , (64a) 
^2g.eo(r.+r„)-i . (64b) 



[^0 - TeFo] 



£0(Fe+F„) 



It is obvious that Eq. H64a|l has the solution ig = 

±-\/f ef o, among which only the positive root satisfies 
Eq. (j64b|l . Substituting this result in Eq. I|64b|) . we ar- 
rive after some algebra at 

\h\2 - -P 1 (l-/3)~ (.rgo/2)(r.+rj-^ /(.t-x 



In the deep Kondo limit, where N ^ — \b\'^ ^ N ^ and Using the previous result, we may define the Kondo tem- 
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perature for each channel as 



41 



K 



T°K 



(1 + /?)^ 
(1-/3)'" 



j^g(7r£o/2)(r,+r„)- 



and obtain the renormalized level position: 



where we are back to the original basis and the trace 
includes also the sum over spin indices. Here, G< is the 
lesser matrix orbital Green function, which is related to 
the advanced and retarded G"" matrix Green functions 
through the expression 



(69) 



V2 



(67) 



Equations (|56|l and H67|) . which are the main results of 
this section, give the Kondo temperatures and level posi- 
tion for arbitrary mixing (3. Note that T^ + To = 2r does 
not depend on /? and therefore the Kondo temperature 
depends on the orbital mixing only through the prefac- 
tor. While Tf^ changes very little with /3, reduces 
down to zero as /3 — > 1 (maximum mixing). Similarly, 
goes from Tk = NT\b\^{(i = 0) = (L>/\/2) exp(7reo/4r) 
to zero, in agreement with the Friedel sum rule. From 
the above results, we conclude that the odd orbital be- 
comes decoupled at maximum mixing, where wc are left 
with SU(2) Kondo physics arising from spin fluctuations 
in the even orbital channel. This SU(4)-to-SU(2) tran- 
sition as mixing increases is illustrated in Fig. 2, where 
the SBMFT parameters are plotted versus (3. 

Now we are in position to calculate transport prop- 
erties. For this purpose, it is more convenient to write 
SBMFT equations in the matrix form: 



1 



^TrG<(.) = l 



(68a) 



A|^P + ^ / ;^Tr{E'-G<(e) + S<G'^(6)} = 0, (68b) 
N / 2m 



with being the lesser matrix self-energy. The explicit 
expressions for these matrices are 



G^e) 



1 



(e - £0 - 
e - £o - iTji 



(70) 



K 

i(3Tk 



and G"' given by direct complex conjugation of G°). The 
lesser matrix self-energy reads 



f Tk I3Tk 



S]< - -2z [hie] + fiiie)] x 
whereas in the same way the retarded matrix self-energy 



(71) 



\PTk Tk 



(72) 



Inserting Eqs. {TUJ and ifTTl) in Eq. (|^ . we obtain the 
lesser orbital Green function: 



G^ 



-iTk 



(e - eo)4 + 2(1 + /32)r|.(e - e^f + {[3^ - 1^ 



(e-eo)'+ 7^(1-/32) /3 [(e - eo)' 
p [{e - So)' ~ Tk{1 ~ /?')] (e - eo)' + Tk{1 - /?') 



-^(1-/32) F-^> 



Using the explicit expressions of the self-energies and the 
nonequilibrium Green functions, we write Eq. I)68bl) in 
the simple form: 



(74) 



Solving in a self-consistently way Eqs. (|68a|l and 174|) 
for each dc bias Vdc, one gets the behavior of the two 
renormalized parameters in nonequilibrium conditions. '^^ 

The electrical current has in appearance the same form 
as the conventional Landauer-Biittiker current expression 



for noninteracting electrons: 



2e f de 
2^ 



r(e,ydc)[/L(e)-/fl(e)] 



(75) 



Here caution is needed for a correct interpretation of 
Eq. (|75|l . since it contains "many-body" effects via the 
renormalized parameters that have to be determined for 
each Vdc in a self-consistent way. As a result, the trans- 
mission T(e, Vdc) possesses, in contrast with the noninter- 
acting case, nontrivial behavior with voltage. The non- 
linear conductance is calculated by direct differentiation 
of Eq. ((75|) with respect to the bias voltage: Q = dl/dVdc- 
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In the limit Vdc — (at equilibrium), the linear conduc- 
tance Qq is given by the well-known expression: 



2e2 

6»^-T,0), 



(76) 



where the transmission is 
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FIG. 13: (Color online) Transition from SU(4) to SU(2) 
Kondo physics as obtained from SBMFT: As the orbital mix- 
ing is increased, the SU(4) Kondo effect reduces to the SU(2) 
spin Kondo effect. This is reflected by (a) the position of the 
Kondo resonance as well as by (b) the reduction of the odd 
Kondo temperature down to zero. See the main text. 



T(e) = Tr{G''f G""f } , 



(77) 



Finally, inserting Eqs. JTDJ) and JT^J) in Eq. |77I), one ar- 
rives at the explicit formula for the transmission 



J 



(e - io)* + 2(1 + /32)r|. (e - s,y + [P^ - l^T^ 



(78) 



r 



which is the main result of this part. It is remarkable 
that the linear conductance Qo does not depend on [3. 
In particular, for the SU(4) Kondo model (/3 = 0), the 
transmission takes the simple form: 



T{e) 



2T2 



(79) 



In this case the resonance is pinned at e = eg = Tk with 
the width given by Tk\ this leads to T(0) = 1 and in 
consequence f/o — 2e^//i. For {3 — 1 corresponding to the 
two-level SU(2) Kondo model, Eq. (|7S)l reduces to 



ATI 



2 ' 



(80) 



which leads to the resonance at e = Eq = and again 
Qa = 2e^//i from T(0) = 1. As pointed out, this fact 
makes both Kondo effects indistinguishable through the 
linear conductance measurement. 

All these features are clearly observed in Fig.^J where 
the transmission for different amounts of mixing, i.e., 
different values of /3 is plotted. For (3 = the trans- 
mission peak is located at Tk as expected, whereas for 
[3 = 1 this moves towards e = 0. During the transition 
from the SU(4) to the two-level SU(2) Kondo model, the 
transmission gets narrower and develops a "cusp" , signal- 
ing competition between even and odd channels. This is 
manifested by the following from of the transmission: 



{l + PfT, 



K 



{l-PfT, 



K 



(e-eo)2 + r|.(l + /3)2 



[e-eoY + Tl{l-(3Y 



(81) 



r 



Note that both channels are resonant at the same energy 
eo but have different widths (Fg and Fq), which explains 
the "cusp" behavior. Here we speculate that finite split- 
ting 5e originating from charge fluctuationa^ (not in- 
cluded at the SBMFT level) would give rise to two split 



resonances for /3 7^ 0, namely, ^ = 5e ± vfVFo- 
This is confirmed in the next section, where we present re- 
sults obtained from full NCA calculations including fluc- 
tuations. Eventually, for (3 = 1 the competition does not 
exist and the transmission does not display the cusp. 
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FIG. 14: (Color online) Equilibrium SBMFT result: Trans- 
mission T(e) as a function of the frequency for several values 
of /3. The left inset displays the Kondo temperature versus /3. 

B. Non-crossing approximation Method 

The SBMFT suffers from two drawbacks: 1) It leads 
always to local Fermi liquid behavior and 2) there arises 
a phase transition (originating from breakdown of the lo- 
cal gauge symmetry of the problem) that separates the 
low-temperature state from the high-temperature local 
moment regime. The latter problem may be corrected 
by including 1/N corrections around the mean- field so- 
lution. The non-crossing approximation (NCA)i2iMi^ is 
the lowest-order self-consistent, fully conserving, and $ 
derivable theory in the Baym sensed which includes such 
corrections. Without entering into details of the theory, 
we just mention that the boson fields in Eq. H54I) , treated 
as averages in the previous subsection {b{t)/^/N — > 
(b) I \/N = b) , are now treated as fluctuating quantum 
objects. In particular, one has to derive self-consistent 
equations of motion for the time-ordered double-time 
Green function (with subindexes omitted): 

zG(i,t') ^ {TJ{t)p{t')), 

iB{t,t') = {TMt)b\t')), (82) 

or in terms of their analytic pieces, 

iG{t,t') = G>{t,t')0{t-t')-G<{t,t')0{t' ~t), 
iB{t,t') = B>{t,t')e{t-t')+B<{t,t')9{t' -t).{83) 

A rigorous and well-established way to derive these 
equations of motion was introduced,^^ and related to 
other non-equilibrium methods such as the Keldysh 
method. Here, we just present numerical results of the 
NCA equations for our problem and refer interested read- 
ers to Refs. I3al43l44ll45l for details. 

In particular, the DOS is given by 

Pie) = -- E MAlAe)], (84) 

v—e.OjfT 
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FIG. 15: (Color online) NCA results: Density of states around 
£ = at T = 0.25Tk (left) and T ^ Tk (right) for several 
values of /3. The inset shows a close-up of the /3 — 0.5 curve 
(thick dashed), together with the individual even and odd 
channel contributions (thin dashed). 

where A^^{e) is the Fourier transform of the retarded 
Greens function: 

a: At) = G:jt)B<{-t) - G<At)B'^{~t). (85) 

The DOS for several values of /? at two different temper- 
atures is plotted in Fig. ^1 Interestingly, the cusp be- 
havior of Fig. ll4l in the previous subsection becomes split 
for the even and odd channels. This is illustrated in the 
inset, where the curve corresponding to /? = 0.5 is plot- 
ted together with the individual even and odd channel 
contributions. As we anticipate, the presence of charge 
fluctuations induces splitting of Sq due to the different 
renormalization arising from different couplings for the 
even and odd channels Fe/o [see Eqs. and 

V. CONCLUSION 

We have considered the single-electron transistor 
(SET) device with the CNT QD or VQD as the cen- 
tral island in the Kondo regime. In particular we have 
examined the case where the CNT QD or VQD has a 
high symmetry so that the orbital quantum numbers are 
conserved through the system. Emphasis has been paid 
on how different Kondo physics, the SU(4) Kondo effect 
or the two-level SU(2) Kondo effect, emerges depending 
on the extent to which the symmetry is broken in real- 
istic situations. Employed are four different theoretical 
approaches: the scaling theory, the NRG method, the 
SBMFT, and the NCA method to address both the lin- 
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ear and non-linear conductance for given external mag- 
netic fields. Our results show that there is no way to 
distinguish experimentally the SU(4) Kondo effect and 
the two-level SU(2) effect by means of the linear conduc- 
tance (which is a low-energy property) alone. The SU(4) 
Kondo physics, which arises with higher symmetry, can 
be observed exclusively only by the non-linear conduc- 
tance (a higher-energy property) in the presence of finite 
magnetic fields, as in the recent experiments^ The sym- 
metry breaking (either the orbital anisotropy 1— r2/ri or 
the orbital mixing Fx) drives the system from the SU(4) 
Kondo fixed point to the SU(2) Kondo fixed point. At fi- 
nite yet sufficiently small symmetry breaking, the SU(4) 
Kondo physics governs the transport in the system at 
relatively high energies (of the order of the Kondo tem- 
perature) while the two-level SU(2) Kondo effect takes it 
over at extremely low energies. This gives another reason 
why the indication of the SU(4) Kondo physics should be 
investigated by means of the non-linear conductance in 
the presence of external magnetic fields. 



with the initial Hamiltonian 



Hn 



/j.m (7 



(A2) 

where the fermion operators ffj,^N,a- have been introduced 
as a result of the logarithmic discretization and the ac- 
companying canonical transformation, A is the logarith- 
mic discretization parameter (taken to be A = 2), 



1- A 



-N 



- A-(2^-l)][l - A-(2^+l)] ' 



(A3) 



and 



Hd = C 



D 



(A4) 
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APPENDIX A: NUMERICAL 
RENORMALIZATION GROUP 

The Hamiltonian allows both charge fluctuations 
and spin fluctuations. Gharge fluctuations (accompanied 
by the particle-hole excitations) occurs at high energies 
while spin fluctuations prevail at low energies. Therefore 
in order to understand the low-energy properties of the 
system, it is useful to take the RG approach and to obtain 
an effective Hamiltonian allowing only the spin fluctua- 
tions. One may follow the three-state perturbative RG 
procedure (scaling theory): One first renormalizes the 
Anderson-type Hamiltonian Q until the charge fluctu- 
ations are completely suppressedii^ (see also Ref. |23) , 
performs the SW transformation^ to obtain a Kondo- 
type Hamiltonian where spin fluctuations are described 
by the spin operators, and renormalizes further the re- 
sulting Kondo- type Hamiltonian. '^^ The RG equations for 
the coupling constants in the Hamiltonian allow one to 
identify physically interesting flxed points and associated 
scaling properties. 

We follow the standard procedure to implement the 
NRG calculations (i2i22i2i and evaluate various physical 
quantities from the recursion relation {N > 0) 

Hn+1 = \^Hn 

+ ^N+l (4,iV,./M,N+l,- + h.c}j (Al) 



with ( — 2(1 + 1/A) ^. The couphng constants V^^.m 
have been defined to be 



fi^ni I 



ttD 



(A5) 



where po is the density of states of the leads at the 
Fermi energy. The Hamiltonian in Eq. (|A1|I has been 
rescaled for numerical accuracy, and the original Hamil- 
tonian is recovered by 



— = lim — — 



(A6) 



with = CA(^"^)/2. At each iteration of the NRG pro- 
cedure, we calculate the local spectral density;^ which 
determines the transport properties through the dot: 

Ara^^, iu) = A>^, (u) - A<^, (w) (A7) 

with 

(A8a) 

ALn':AE) = - E "^I'a !«) («l d,na |0) 5(i? + - Eo) 

(A8b) 

where \a) represents the many-body state of the system 
with energy Ea (with a = corresponding to the ground 
state). 
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